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NODAL SETS OF SCHRODINGER EIGENFUNCTIONS IN 
FORBIDDEN REGIONS 

YAIZA CANZANI AND JOHN A. TOTH 


Abstract. This note concerns the nodal sets of eigenfunctions of semiclassical 
Schrodinger operators acting on compact, smooth, Riemannian manifolds, with no 
boundary. We prove that if H is a separating hypersurface that lies inside the 
classically forbidden region, then H cannot persist as a component of the zero set 
of infinitely many eigenfunctions. In addition, on real anal3rtic surfaces, we obtain 
sharp upper bounds for the number of intersections of the zero sets of the Schrodinger 
eigenfunctions with a fixed curve that lies inside the classically forbidden region. 


Let {M,g) be a smooth, compact, Riemannian manifold with no boundary. Write 
Ag for the Laplace operator, and given any smooth potential V G M), consider 

the Schrodinger operator acting on Lp‘{M) dehned as 

P{h) = -h^Ag + V, 

where h G (0,1]. Let E G M be a regular value for the total energy function p{x,^) = 
defined on T*M, and write for the the classically forbidden region 

Ob := {x G M : V{x) > E}. 

In this paper we study the nodal sets of Schrodinger eigenfunctions (with energy close 
to E) inside the classically forbidden region, in the semiclassical limit /i —)• O"*". Consider 
L^-normalized Schrodinger eigenfunctions with 

P{h)(f)h = E{h)(j)h and E{h) = E + o{l) as h —)• 0"^. (1) 


There is a large literature devoted to the study of the zero sets of Laplace eigen¬ 
functions, 

= {x € M : <f)h{x) = 0}, 

on compact manifolds. We refer the reader to |Z2j for a detailed list of references. 
The Hausdorff measure of the zero sets, their distribution properties, the number 
of nodal domains and their inner radius, have been extensively studied (although 
many open problems remain, even for surfaces). More generally, it is natural to study 
the properties of zero sets of Schrodinger eigenfunctions inside the classically allowed 
region where V < E. Many of the known results in the homogeneous case where 
V = 0 extend to Schrodinger eigenfunctions in the allowable region (see [.Tin| ). In 
contrast, very little is known about the zero sets of Schrodinger eigenfunctions inside 
the classically forbidden region where V > E. In dimension one, it is known that the 
eigenfunctions of the Harmonic Oscillator have no zeros in the forbidden region and 
in recent work, Hanin-Zelditch-Zhou [HZZ| have proved that in any higher dimension 
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the expected value of the measure of the zero set of random eigenfunctions of the 
harmonic oscillator inside any ball is of order . We are not aware of any other 
results addressing the behavior of zero sets of Schrodinger eigenfunctions inside the 
classically forbidden region. 

Our first result addresses the issue of nodal persistence: Can a fixed hypersurface 
H be contained in the nodal set of an infinite subsequence of eigenfunctions? This 
question was answered on the flat torus by J. Bourgain and Z. Rudnick in [BRj . 
They proved that if R = 0 and H C T” is a hypersurface with non-zero principal 
curvatures, then H cannot lie within the zero set of infinitely many eigenfunctions. 
On general manifolds with V 7 ^ 0, we prove that no embedded separating hypersurface 
contained entirely within the forbidden region can persist as part of the zero set 
for infinitely many eigenfunctions. 




Figure 1. Level sets of eigenfunctions of —h'^Ag -|- R on a square torus, 
where E = 1 and R is a periodized sum of two bumps 4e“i°(R+o.3)^+(y+o.3)^) _|_ 
3g-i5((a:-o.6) +{y-0.7) pictures Correspond to h = 0.01 and h = 0.005 

respectively. Tones of blue describe negative values, tones of red describe 
positive values. In the plot the value zero (and very small values as well) are 
depicted as white. 


Theorem 1. Let {M,g) be a smooth, compact, Riemannian manifold with no boundary 
and let V G C°°{M). Consider a sequence {cfh} as in ([^. Suppose that H C Qe is 
an embedded separating hypersurface that encloses a bounded domain contained in Qe- 
Then, there exist constants Ch > 0 and ho > 0 such that 

\\4>h\\L^H) > and \\du4>h\\L^H) > 

for all h G (0, Lq], In particular, for every subsequence {hj}j with hj —)• 0 as j —)■ 00 , 
there is an integer jo > 1 with the property that for all j > jo 
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Remark 1. Theorem extends to the case where M = provided that V G 
extends holomorphically to a complex wedge domain = {z € : 

|Imz| < ^(z)} and that it satisfies V[x) > C'{x)^ for some k G Z+ as |x| —)• oo. Here, 
C and C are positive constants. 

Assume from now on that (M, g) is a compact, real analytic surface and let R <Z^e 
be a real analytic closed curve that bounds a region inside Unique continuation 
results like the one in Theorem have important implications for the study of asymp¬ 
totic oscillation properties of eigenfunctions, including estimates for the intersection 
number (1 H} (see for example |TZ[ IHTl ICTj i. Our second result is a deter¬ 

ministic upper bound for the nodal intersection with a fixed real-analytic curve H on 
a Riemannian surface contained in the classically forbidden region. 


Theorem 2. Let {M,g) be a eompact, real analytic surface with no boundary. Let {(f^h} 
be real valued eigenfunctions satisfying 0, where we also assume that the potential 
V is real-analytic. Suppose that H <Z VLe is a simple, closed, real analytic curve that 
bounds a region inside Qe- Then, there exists Ch > 0 and /iq > 0 such that 

for all h G (0, /iq] . 

To prove Theoremj^we use the restriction lower bound in Theorem [^together with a 
potential layer formula for the eigenfunctions inside the forbidden region. Bounding the 
number of zeros on the fixed curve is then reduced to estimating the complexihcation of 
a particular Green’s operator in the forbidden region. We control the complexihcation 
of the Green’s operator using off-diagonal decay estimates for the real kernel (see 
Proposition]^ together with h-analytic Gauchy estimates recently proved by L. Jin in 

g^. ^ 

For individual eigenfunctions, one can see that the CHh~^ bound in Theorem]^ is 


sharp on surfaces of revolution (see Section 2.3) and agrees with the upper bound in 
Yau’s conjecture [Y1,Y2] for nodal volume in the homogeneous case. Nevertheless, 
it is reasonable to expect that in many cases one should be able to improve on this 
bound. As mentioned, for random eigenfunctions of isotropic harmonic oscillators, 
the computations in jHZZ] show that the expected value of the nodal lengths in the 
classically forbidden region are of order Gonsequently, at least for random 

waves, it is reasonable to expect generic intersection bounds of the form = 

in the case for which the forbidden region is unbounded. We hope to return 
to this question elsewhere. 

0.1. Organization of the paper. In Section we prove Theorem using an el¬ 
ementary argument with Green’s formula and quantitative unique continuation for 
the eigenfunctions. In Section we study nodal intersection bounds by reproducing 
the eigenfunctions in the forbidden region using a suitable Green’s operator whose 
complexihcation we need to control. Assuming that we have suitable bounds on the 
complexihcation of the Green’s operator we then prove Theorem [fusing the restriction 
lower bound in Theorem In Section 2.3, we show that the upper bound in Theo¬ 
rem is sharp. In Section we give a detailed analysis of the kernel of the Green’s 
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operator on a compact manifold. In particular, we show that the kernel can be locally 
complexified away from the real diagonal {(x,x) G M x M} and obtain exponential 
decay estimates in h for the complexified kernel. 

0.2. Acknowledgements. The authors would like to thank Christopher Wong for 
sharing his MATLAB code to compute estimations of Schrodinger eigenfunctions on 
the square torus. 


1. L^-lower restriction bounds 

We note that because of the quantum tunnelling effect, the wave functions are known 
to have positive mass inside the classically forbidden region. Indeed, by Carleman type 
estimates |Zwl Theorem 7.7], for every open set U <Z^e there exists a positive constant 
C = C{U) > 0 for which 

> e as /i—(2 ) 

The result in Theorem is an analogue of the lower bound in ([^ for the eigenfunction 
restricted to a hypersurface H G and is a crucial step in the proof of Theorem 
The only condition that we impose on H is that it must bound a domain that is 
contained entirely inside ^Ie- 

We note that the exponential lower bound in Theorem is quite delicate since 
despite the fact that inside the forbidden region the eigenfunctions have positive mass, 
they are exponentially small in h. Indeed, consider the Agmon metric Qe = (V —1^)+9 
and associated distance function dE{x) := dist g^ {x, M\Qe)- By the standard Agmon 
estimates |Hell Proposition 3.3.4], it follows that for any e > 0, 

\dx4>h{x)\ = Oe.a (e "h ^ (3) 

locally uniformly in x G fie. In particular, given a smooth hypersurface H C ^e, h 
follows from Q that for dE{H) := min{d^(x) : x G M}, one has 

\\(l^h\\L'^{H) = 0{e h y (4) 

We may then view Theorem as a partial converse to @ under the assumption that 
H C He is a separating hypersurface. 

l.I. Proof of Theorem Let H C i^E be a separating hypersurface that bounds a 
smooth domain Mh C He. Since Mh C fls and E{h) = E + o(l) as /i —>• 0, it follows 
that if /lE > 0 is sufficiently small, then there exists Ce > 0 so that 

V{x) — E{h) > Ce for all x G Mh (5) 

and all h G (0, He]- 
By Green’s Theorem, 

I dVg T / (/l /\g)(p}i (f)}idVg - h I di/ijifi (pjidCTg, 

J Mjj J J dM^ 

where v is the outward normal vector and ag is the induced volume measure on dMn- 
Thus, since —h'^Agcjih + Vcph = E{h)(j)h, it follows that 
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g4>h\\‘L‘2{M„) + ((^ “ E{h))(j)h, 4>h)L2{MH) — h‘^{du4>h, (t’h)L^{dM„)- 

Using the non-negativity of and ([S]) we obtain that for all h G (0, He] 

<h?{du(t)hAh)L'2(dMH)- (d) 

An application of the Cauchy-Schwarz inequality in Q gives 

CEh~‘^\\(t)h\\\2(^MH) - \\^h\\L^{H) \\du(l)h\\L-2(H)- ( 7 ) 

By the unique continuation lower bound Q there exists C > 0 with \\4>h\\\2(^]\^^^ > 
g-C”//! £qj. ^ small enough and therefore there exists He > 0 such that 

CEh~‘^e~h < \\cl)h\\L'^(H)\\dv(t>h\\L'^(H), 

for all h G (0, /ig]. Theorem then follows from the Agmon estimates in Q. 

□ 


2. Nodal intersection bounds 


Here we present the proof of Theorem]^ (Section 2.1), and show that the upper 
bound on the number of nodal intersections is saturated for surfaces of revolution 
(Section |2.3[ ). 


2.1. Proof of Theorem We continue to assume that C but here we make 
the additional assumption that the Riemannian manifold (M,g), the potential V and 
the hypersurface H are all real-analytic. Also, in the following, we are only interested 
in the case in which M is a surface. Let q : [0,27r] —)• iL be a C‘^, 27r-periodic, 
parametrization of H. To bound the number of zeros oi 4‘h o q : [0, 27r] —)• M we 
consider its holomorphic extension ((ph o q)^ : —?■ C to the complex strip 

= {t G C : Ret G [0,27r], |Imt| < r} 

for some r > 0, and use that n H} < H H^}. Then, the zeros of 

{4>h o q)'^ are studied using the Poincare-Lelong formula; 

ddlog\{(t)hoq)'^{z)\‘^ = 

According to m Proposition 10], there exists C > 0 so that 

#{Z^,nF} <Cmaxlog|Ff(t)|, (8) 


where (t) with t G is the holomorphic continuation of the normalized eigenfunc¬ 
tion traces 


p _ <t>h{q(t)) 
UhWL^iH)' 


(9) 


Note that by Theorem we know that ||(/>h||H > e £qj, ^ ^ (0,/io] with Hq 

sufficiently small, and this implies that Th(t) is well defined. 
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It follows that we shall need to control the complexification F^(t) to obtain upper 
bounds on n H}. Without loss of generality we assume that H C 

where Qh C fie is a domain containing 7 , a closed C‘^ curve, as its boundary. Our 
goal is to find a double layer jumps formula that reproduces (phix) for x € H m terms 
of its values along 7 . 

Let X G C^{M, [0,1]) with xix) = 1 for all x G Qh and with suppx C Qe- Consider 
the auxiliary global metric given by 


+ xGM. (10) 

From now on, to simplify notation, we simply write g^^ for ^ since the dependence 
of the latter on h is only in the constant eigenvalue term and is of no real consequence 
as far as h-pseudodifferential calculus is concerned. 

Then, since n = 2 , it follows that = [(C — E{h))x + (1 — x)]~^ and so, in 

particular, for all x G 

(-h^^gn +'^)4>h{x) =-h^{V{x) - E{h))~^Ag4>h{x)+ (l)h{x) = 0. (11) 

iiE 

We consider the Green’s operator 


G{h) = + 1) 


-1 


Since {—h^Ag^ + l)G{x,y, h) = 6x{y), (11) implies that for all x G Mh 

4>hix) = 



,y 


+ '^)G{x,y,h)(j)h{y)dv{y) 



Gix,y,h){-h^Ag^ 


+ ^)4>h{y)dv{y) 

( 12 ) 


By Green’s formula, it then follows that for x G 12/^, 

(phix) = G(x, y, h)dyy(t)h{y)da{y) - J duyG{x, y, h)(j)hiy)da{y)^ , (13) 

where Vy is the outward normal vector at y G 7 . Let r : [0,27r] —?■ 7 be a C"^ 
parametrization of 7 . Restriction of the outgoing variable a: in (13) to LI yields the 
potential layer formula 

4>hiq{t)) = h'^ ( G{q{t),r{s),h)d^y(ph{r{s))da{s)-h‘^ f duyG{q{t),r{s), h)4>h{r{s))da{s) 
J 'y J 'y 

(14) 

In order to control Ef^{t) in (|^, one needs an upper bound for the holomorphic 
continuation {(f>h ° q)'^- The latter amounts to estimating the complexihcation of (14). 

Given r > 0, let q'^{t) denote the holomorphic continuation of the parametriza¬ 
tion y : [0, 27r] —)■ LL to the strip We claim the following result. 


Proposition 3. Suppose dg{H,Hs) > e for some e > 0. Then, there exist constants 
C{e) > 0, r(e) > 0, ho{e), such that for h G (0,/ 10 (e)], 

\G^{q^{t),r{s),h)\=0{e~^) and \dyyG^{q^{t),r{s),h)\ = 0{e~^), 

(15) 

uniformly for (t, s) G x [ 0 , 27 r]. 
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Proposition!^ is a consequence of a more general result, Theorem]^ which we prove 
in Section (see Remark [^. In Theorem]^ we show that the kernel of the Green’s 
operator G{h) along with its derivatives can be locally complexified off-diagonal main¬ 
taining the exponential decay exhibited in the real domain. 

Substitution of the estimates in Proposition into the complexification of (14), 
combined with an application of Theorem and the Cauchy-Schwarz inequality gives 
the existence of positive constants C, ho and dn such that 

\F^it)\ < (16) 

\\\(Ph\\L2{H) \\(Ph\\L2{H) J 

for all h G (0, ho]- Then, by ([^, there exists Ch > 0 such that 

#{Z^^r^H]<CHh-\ (17) 

as desired. Assuming that Proposition holds, this concludes the proof of Theorem 

E □ 

2.2. Estimates in classically allowable versus forbidden regions. It is inter¬ 


esting to contrast the growth estimates in (16) with the case where H is contained 
in the classically allowable region. For example, when H C ^ where 11 C is a 
piecewise-analytic planar domain, and cph is a homogeneous eigenfunction (satisfying 
either Neumann or Dirichet boundary conditions), one can show that (see m Lemma 

11 ), 

I'f (t)| < 

y \\(Ph\\L'^{H) \\(Ph\\L'^{H) J 


\F^ 


(18) 


In (18), the constant C 2 = where dP is the complex- 

ifieodistance between dll and H, and q, r are as defined in Section |2.1[ Thus, C 2 


is positive in contrast with the negative constant —Ci{e) appearing in the forbidden 
case. This is due to the fact that the Green’s kernel G{q{t),r{s), h) is a semiclassi- 
cal pseudodifferential operator that decays exponentially off the diagonal and so does 
the corresponding local complexification G'^{q'^{t),r{s).,h) (see Theorem]^. In the 
allowable region, the Green’s kernel G{q{t),r{s), h) is replaced with the restriction of 
the free Helmholtz Green’s kernel G-^ 2 {q{t),r{s), h) in which has the WKB asymp¬ 
totics G]k 2 (g(l), r(s), h) {2TThyPe’‘PP^')’'"P'iP^{ao{t, s)-|-ai(l, s)h+- ■ ■), provided 

hF(^q r^^jjy^QQd{q,r) > 0. This is the kernel of an /i-Fourier integral operator and the 
phase factor blows up exponentially in /i —)• O'*' upon complexification in q{t), 

unlike in the /i-pseudodifferential case where there is off-diagonal exponential decay in 
h. In view of the Jensen-type growth estimate in ([^, it follows that the constant in the 
OH{h~^) intersection bound in the forbidden region is smaller than the one for the al¬ 
lowable region. However, as the next example shows, the h~^-rate cannot be improved 
in general. We hope to return to discuss these issues in more detail elsewhere. 


2.3. The example of a convex surface of revolution. Here we show that the 
upper bound in Theorem[^is sharp. To do this, consider a convex surface of revolution 
generated by rotating a curve y = f{r) about the r-axis with / G C'"'([—1,1],M), 
/(I) = /(—I) = 0 and in addition require that f''{r) < 0 for all r G [—1,1] so 
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that the surface is strictly convex. Let M be the corresponding surface of revolution 
parametrized by 

[-1,1] X [0,27r) 

(r, 0) i-A (r , f{r) cos{9), f{r) sin(0)). 

In these coordinates M inherits a Riemannian metric g given by 

g = w‘^{r)dr^ + f‘^{r)d0‘^ 

where we have set w{r) := Consider on {M,g) an analytic potential 

V (r, 6) = V(r) independent of the angular variable, and let E be a regular energy level 
for the hamiltonian corresponding to —h'^Ag + V with mini/ < E < maxi/. We shall 
construct a curve H contained in the forbidden region {V > E} and a sequence of real 
valued solutions {4’h^.}k of {—h?Ag + V)(j)h^, = E{h)(j)hk where E{hk) = E + o(l) as 
k —)• + 00 , so that 

Consequently, the 0{h~^) bound in Theorem]^ is sharp. 







Figure 2. Illustration of the level sets of an eigenfunction in an el¬ 
lipsoid where the potential is concentrated close to the poles. The 
colouring scheme is the same as in Figure 1. 

We seek eigenfunctions of the form 9) = Vh{r)ijjh{9) that solve 

{-h^Ag + V)(l)h = E{h)<t>h. 

Since the Laplace operator in the coordinates (r, 9) takes the form 


1 


d 


A - , fjr) d 

^ w{r)f{r) dr \w{r) dr 

we have that the functions Vh and must satisfy 




w‘^{r) d^ 
P{r) d9‘^ ’ 


d^ 


- ^'^h{9) = mli^h{9) 


(19) 


and 
- h 


2 fir) d ( f{r) d 


rhir) +/ (r)(l/(r) - E{h))vhir) = -mih^w^{r)vhir), (20) 


w{r) dr \w{r) dr 

for some nih G Z. From now on let {/ifcjfc be a decreasing sequence with hk 


O"'' as 

-|-oo and such that mh,, = l//ifc G Z. One can choose the solution of (19) to be 

= (21) 
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To deal with the radial part one makes the change of variables s i—?• r{s) = 
/o^ ™ follows that 

£ 


-hi 


ds"^ 


+ f iris)){V{r{s)) - E{h)) + w\r{s)) Vh^{r{s)) = 0. 


( 22 ) 


To finish the argument, let ro G {r G [—1,1] : V{r) > E} and set HrQ to be the 
curve 

Hro := {ir,e) : r = ro}. 

Since cos{9/hk) has 2/i^^ zeros for 6 G [0, 27r), our claim is established once we set 


9) 


with Vhv. a solution to (22). 


Vh,{r)) 


3. The Green’s operator G{h) and its complexification 

Let {M,g) be a compact, real analytic Riemannian manifold of dimension n. We 
consider here the associated Green’s operator 

G{h) = {-£Ag + l)-i : C“(M) ^ C“(M). 

The purpose of this section is to study the complexification of the kernel G{x, y, h) 
in the outgoing variable x. Before we state our main result (Theorem]^, we briefly 
review some of the complex analytic geometry that is needed in the formulation and 
proof of Theorem 

3.1. Grauert tube complexification of M. By a theorem of Bruhat-Whitney, M 
has a unique complexification with M C M'^ totally real that generalizes the 
complexification of ME to C”". The open Grauert tube of radius e is defined to be 

= {z€M^: ^g{z) < r}, 

where on is the unique solution to the complex Monge-Ampere equation. 
For example, in the simplest model case when M = M? and = C^, one has 
^/p^{z) = 2|Imz|. There is a maximal Tmax > 0 for which is defined [LGS) Thm 

1.5], and is a strictly pseudoconvex domain in for all r < Tmax- 

For all r < Tmaxj we identify the radius r ball bundle {BM)t- C TM with {B*M)r C 
T*M using the Riemannian metric. For x G M and 0 < r < inj(M, gf), we let exp^. : 
Bx{0, r) —)• M be the geodesic exponential map defined on the geodesic ball ^^(O, r) C 
T*M. We denote the lifted exponential map to all of {B*M)t by 

Exp : {B*M)r —)• M, Exp(x,.^) = exp 3 ,(^). 

Since {M,g) is real-analytic, for hxed x £ M and 0 < r < injM, the geodesic exponen¬ 
tial map exp^, : Bx{0, r) ^ M admits a holomorphic continuation exp^ ; {Bx{0, r))*^ —>• 
M'^ in the fiber ^-variables. For 0 < r < Tmax, we define the associated complexified 
lifted map by 


Exp^ : {B*M)r M^, Exp^(x,0 = exp^(fO- 
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The map Exp^ gives a diffeomorphism between {B*M)r and with the property 
that (Exp'^)*(/9g) = I 'Ig. Consequently, {B*M)r = as complex manifolds via Exp®'. 
Also, the map 

7rJ^^(Exp®^(a:,0) = x (23) 

is an analytic fibration. The hbers correspond to imaginary directions over 

the totally real submanifold M C 

In general, we shall denote by (jP : —>■ C the complexification of a real analytic 

function (/> : M —>■ C. 

Eix xq G M. The map 

Bxo (0, r) M 

V = r{x) eA exp2,^(r/) = x, 

is real analytic near the origin and so it can be holomorphically extended as 

{BxP0,r)f ^ ^ {B*MP 

V + iC = f{x,0 ^ exp^^(r/ + 0 = (x,0- 

By Lemma 1.18 in [LGSj . this coordinate system satishes f{x, 0) = r{x) and /(xq, 0 = 
Identifying the point (x,^) G B*M with exp^(i^) G as described above, one 
has TTj^Px,^) = (Exp®'(x,^)) = x = exp^^{r]). In view of Lemma 1.18 of [LGSj we 

will use holomorphic coordinates (t/, C) = /(x,0 the complex manifold {B*M)t-. 

Erom now on, in a coordinate neighbourhood of xq G M, we write z = z{x,^) G 
for complex coordinates where 

z = (Re z, Im z) 

with 

Rez := Re/(x,^), and Imz := Im/(x,^). (24) 

Note that with this notation t^;^{z) = Tr^j{x,^) = exp^^{r]) = exp^^(Re f{x,^)) = 
expj.^^ (Re z), and so ttm{z) is identified with Re^;. 

We also have that G C°°{M^) is a strictly plurisubharmonic exhaustion function, 
and by Taylor expansion around Imz = 0 it follows that p‘^{z) = 4|Imzp + 0(|Im2;p). 
Thus, for 2 : G M!^ the function ^J~p^{z) ~ 2|Im2:| and it will sometimes be convenient 
to work on the subdomain of given by 

G : |Im 2 ;| < 

with C > 0 sufficiently large. 

3.2. Statement of the main result. Let e > 0 and consider the e-diagonal neigh¬ 
borhood A(e) = {(x,y) G M X M : dg{x,y) < e}. Under the analyticity assumption 
on {M,g), 

G{-, •, h) G C^^iM X M\A(e)). (25) 

Indeed, for any h G (0,1), —h?/S.g -|- 1 is a real-analytic partial differential operator 
that is uniformly elliptic in /i G (0,1) and for hxed y ^ M and any e > 0, 

{-h'^Ag^x + l)G'(x, y, h) = h~^6{x - y) = 0 


( 26 ) 
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whenever {x,y) G (M x M)\A(e). Consequently, (25) follows by analytic hypoellip- 
ticity. We may then consider the complexification y,h) in the outgoing variable 

for {z,y) G ^T(e) ^ ^ with dg{'K^jZ^y) > e, where r(e) > 0 is sufficiently small de¬ 
pending only on e. The purpose of this section is to prove the following asymptotic (in 
h) supremum bound for the holomorphic continuation G‘^{z,y,h) in a small complex 
Grauert tube over the real off-diagonal domain. 


Theorem 4. Given e > 0 there exists a constant r(e) > 0 such that the Greens kernel 
admits a holomorphic extension G'^{z^ y, h) to (z, y) G x M with dg{'K^jZ, y) > e. 

Moreover, for any fixed Nq G N”, and a G N"" with |a| < Nq, there exists G = 
G{e, Nq) > 0 such that as h ^ 


\dyG^{z,y,h)\ = 0{e , 


(27) 


uniformly for {z,y) G x M with dg{'Kj^z,y) > e. 

We prove Theorem]^ in Section [3.3[ 

Remark 2 (Proof of Proposition]^. Propositionfollows directly from Theorem]^ 
as a special case. Indeed, since dg{H,'y) > e for some e > 0, by choosing the Grauert 
tube radius r(e) > 0 small enough, the estimate in ^7\ i is satished by the Green’s 
kernel associated with the extended Agmon metric defined in (10). We note that 
although is only globally G°° on M, it is real-analytic in the forbidden region ^Ie- 
Thus, the Schwartz kernel, G{x,y,h), of G{h) = {—hfNg^^ -|- 1 )“^ is real-analytic for 
{x,y) G Ob X Q.E- 


To prove Theorem we must first describe the real kernel G{x, y, h) in detail. 
In particular, we prove exponential off-diagonal decay estimates for the real kernel 
G{x, y, h) in Proposition below. Although such results are known to experts, we 
could not find a reference in the literature containing all the details we need here. As a 
result, for completeness, we carry out in detail the /i-analytic parametrix construction 
in using the method of analytic stationary phase, keeping careful track of the 
various remainder terms. 

Consider the cut-off function 


X{x,y) ■■= x{dg{x,y)), (28) 

for X G C'^([0,-hoo]; [0,1]), with suppx C [0, inj(M, 5 )] and 

X{t) = 1 for t G [ 0 ,inj(M, 5 r)/ 2 ] . 

We cover the manifold with coordinate patches so that if xix,y) 7 ^ 0 then x and y 
belong to a common coordinate neighborhood. 

Proposition 5. The Schwartz kernel of the Greens operator admits a decomposition 
of the form 

G{x, y, h) = Ag{x, y, h) + Rg{x, y, h). (29) 


Here, 


AG{x,y,h) 


x{x,y) 

(2-71/i)” 
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with aGix,y,r],h) = Ylk=o^ h^Wk{x,y,r]) G ^ for some Cq > 0. Also, for each 
a, /3 G N** there exists Ca^/s > 0 and ho = ho{a, /3) > 0 such that for all h G (0, ho], 

\d^d^RGix,y,h)\=0{e-^-^^/^), (30) 

uniformly in x,y G M. 


Remark 3. We note that writing Gy{x]h) := G{x,y,h) it follows that for x G 
M\By{e) the function Gy satisfies P{h)Gy{x) = {—h'^Ag + l)Gy{x) = 0, and since 
P{h) is uniformly elliptic for h < ho, it follows from the maximum principle that 

max \Gy{x)\ < max |Gy(a:)| = 0{e~‘^/^), 

{x:dg[x,y)>\ inj{M,g)} {x. dg{x,y)=\\n]{M,g )} 

where the last estimate follows directly from Proposition Thus, it suffices to 
bound x{x,y)G{x,y,h) since the far off-diagonal part of the Green’s kernel (1 — 
x{x, y))G{x, y, h) is controlled by the former and is absorbed into the remainder term 
RG{x,y,h) in Proposition]^ 

We prove Proposition in Section 

y, 

Ag{x, y, h) = aG(x, y, rj, h) dy, 

(ivrh)" Jku 

and so, restriction of Ag{x, y, h) to points x,y ^ M with dg{x, y) > s > 0 gives 

/ In 

Re [i{gy^ exp-\x),r]) - j <- —{rj) 

for some C > 0. As a consequence of Proposition one gets the off-diagonal expo¬ 
nential decay estimates for G{x,y,h). 

Corollary 6. Let e > 0 and No G N. Then, there exists G = C{e,No) > 0 such that 
for all a, /3 G N"’ with |a| < No and |/3| < Nq, 

\d^d^G{x,y,h)\ = O (31) 

as h —)• 0^, uniformly for {x,y) G M x M with dg{x, y) > e. 

In order to avoid breaking the exposition at this point, we defer the proof of Propo¬ 
sition to Section |3.4| and proceed with the proof of Theorem 


3.4 In terms of normal coordinates centered at 


3.3. Proof of Theorem For y G M let Gy{x, h) := G{x, y, h) be the real Green’s 
function in (25). Fix e > 0, and consider the ball 

By{e) = {x G M : dg{x,y) < e}. 


Then, since —mAg -|- 1 is h-elliptic, in view of (26), by the semiclassical Cauchy esti¬ 
mates [.Tinl Theorem 2.6], for each xo G M\By{e) there is a coordinate neighborhood 
U C {AI\By{e)) with xo GlA and a positive constant Go, such that for all x G W and 
a G N", 

\dfGy{x, h)\ < G[^\h-^ + |a|)H ||G,( •, h)]]^^^,), (32) 
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for all h G (0,1). Moreover, the estimate (32) is locally uniform in y G M. Let z G 
for a tube radius r > 0 to be determined later. By Taylor expansion around Rez €U, 
and using the Cauchy estimates (32), we have 


\G^{z,h)\< T 

^ a\ 


|l“l 


| a |=0 


<\\Gy{-,h)h^(U)[Y.^0 

\ a \=0 


|a| 


(/i-i + |a|)H 


a\ 


|Imz| 


(33) 


Let Cl > 0 be a large constant to be determined. Splitting the RHS of (33) into two 
terms, we get 

\G^{z,h)\ < Ti{z,y,h) + T 2 {z,y,h), 


for 


Ti{z,y,h) 

T2{z,y,h) 


[( Cih ) ij \ a \, j_i l^ldal 
| a |=0 

\\r^ / i,Mi C*!) + |a|)l“l I I 

\\Gy{-,h)\\L^(^U) 2^ -limbi' 


|a|>L(Ci/i)-lJ 


To control the term Ti in (34), we use the basic estimate 
L(Ci/i)-iJ |aL, _1 I olal 

y ^0 +l«ll |i^ z\H < gCo(l+Cfi)|Im^|//i^ 


E 

| o |=0 


a\ 


(34) 

(35) 


and get that 


|ri(z,y,h)| < ||G,(-,h)|Uoo(^)e^°(i+^")|i”^l/2 


As for the second term T 2 in (35), we use the Stirling-type lower bounds a! > 
(ne)“l"l|a|l"l for all a G N” |.Tin[ (2.9)] to get 

|r2(^,y,h)| < \\Gy{-M\L^(u) E |C2-e|H|imz|H <C3||Gj,(-,h)|Uco(w), 

\ a \> h-'-/Ci 

for some G 3 > 0, provided z G with jlmzj < 

By the real off-diagonal estimates in Corollary]^ there is a constant G(e) >0 such 
that 




(36) 


Since M\By{e) is compact, the theorem follows from (36) and the local bounds for 
|Ti(z,y,/i)| and |T 2 (z,y, h)| above, by choosing jlmzj < r(e) with r(e) > 0 sufficiently 
small. Higher derivatives |(9“G®'(z, y,/i)| are bounded in a similar fashion. □ 
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3.4. Construction of G{x, y, h): Proof of Proposition In this section we prove 
Proposition by constructing an /i-analytic parametrix G{h) for the Green’s operator 
G{h) following closely the treatment in Section 1]. 

Given a compact, real-analytic Riemannian manifold {M,g), consider a complex 
neighborhood of T*M of the form 

{T*M)^ := |(2,C) : ^ G M^, |ImC| < ^(C)} 

with C > 0 fixed sufficiently large. Here, we use the usual convention 

(«C) := + l«sP- 

Following (^, we write a G S^’^{T*M) provided that for all p, g G Z_|_ 

dPdIa = 

uniformly for (x,^) G T*M. We write a G if a ~ h~"^{ao + hai + ...) in the 

standard G°° sense. The symbol a(x, h) is classical analytic (ie. a G pro¬ 

vided a{x,^,h) extends holomorphically to {T*M)'^ and the continuation (denoted by 
a^{x,^, h)) satisfies the following estimates: 


(i) 


-h- 


(O/Coh 

E 

k=0 


kC 




(n) \af{x,0\<CoG^jl{0’^-^ 

(in) (37) 

uniformly for (x,^) G {T*M)'^. In (37) the constant Cq > 0 is sufficiently large and 
Cl > 0 depend on Cq. In the analytic case, we henceforth write a ~ /i“™'(ao + hai -|-...) 
provided ( [37| ) (i)-(iii) hold. 

Gonsider the phase function 


4>{o^x,a^,y) = -{exjp^liy) , ^ (f {ax,y){a^) 

^ z 

for {ax,a^) G T*M and y G M with dg{ax,y) < ini{M,g), and set 


(p*{a, x) = (p{a, x). 


Gonsider the cut-off function 


p{x,y) := p{dg{x,y)), (38) 

where p G C“([0, -|-oo]; [0,1]) is a smooth cut-off function with suppp C [0, inj(M, g)/4\ 
and 

/9(f) = I for f G [0,inj(M,c/)/8] . 

3n n _2 

Given an elliptic symbol b G , consider the corresponding operator Sb{h) : 

G^{T*M) -g G°°{M) 


Sb{h)u{x)= f e p(^ax,x)b{a,x,h)u{a)da. 

Jt*m 


(39) 
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3n n 

Since P{h) = —h'^Ag + 1 is /i-elliptic, there exists q elliptic with 

P{h) O Sb = Sg. 

Using that q is elliptic, one can construct an /i-analytic FBI transform T„(h) : C°°(M) —)• 
C^{T*M) of the form 



Tav{a;h)= / p(ax,y)a{a,y; h)v{y)dvg{y) 

Jm 

(40) 

3n n 

with a G 

satisfying 



Sg{h)oTa{h) = I + Rab{h) 

(41) 

where 

\d^d^Rab{x,y,h)\ = 

(42) 

Consequently, 

P{h)oSb{h)oTaih)=I + Rab{h). 

(43) 


It follows from the parametrix construction ( |43[ ) that the Greens operator G{h) = 
{-h?Ag + l)-i e Oph{S°l~^) is given by 

G{h) = Sb{h) o Ta{h) o (I + Rab{h))-^ 

= G{h) + R{h). (44) 


Here, we have set 

G{h) := Sb{h) o Ta{h) 


and R{h) = -G{h)Rab{h){I + Rab{h)) 


-1 


To compute G{h), we note that from (39) and (40), 

G{x, y] h) = xix, y) f a^,x, y, h) da^da^ 

Jt*m 

where the phase function is 


(45) 

(46) 


ag, X, y) := (/-(aa,, a^,y) - (j)*{ax, a^,x) (47) 

= (exp“J(x) - exp"J(?/) , ^ {dl{ax,x) + dl{ax,y)) 

and the amplitude is 


c{ax, a^, X , y; h) := a{ax, a^,y, h)b{ax, a^, x, h)p{ax,y)p{ax, x). (48) 

We note that the prefactor x{x, y) can be added in ( [i^ since x(x, y) = 1 whenever 

p{ax,y)piax,x) / 0 . 

Given x,y & M with d{x,y) < let a% = Q!^(x,y) be the aa;-critical point 

of the phase <h. We claim that 


al{x,y) = exp 


expyi(x) 


Indeed, in normal coordinates centered at 


dajxn.^{ax,a^,x,y)\^^^^, = -da^ {dl{ax,x) + d^iax^y)) (a^) 


= 0 . 
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Because da^dg{ax,x) = —2exp^^(x) and da^g'^^{a%) = 0, it follows that 

-2{x + y){a^)a- = 0 . 

Therefore, in normal coordinates centered at a%, this gives x = —y. Thus, the critical 
point a% = a^ix, y) is the midpoint of the geodesic segment joining x and y as claimed. 


Since the kernel of G{x,y; h) in (46) involves an integral with amplitude supported 
in the set {ax G M : dg{y, ax) < 1/2}, the analysis is local and from now on we work 
in local coordinates. Then, 


where 


G{x,y,h) = x{x,y) I{a^,x,y,h)da^, 


I{a^, X, y, h) := [ x, y; h) dax 

JR^ 


(49) 


To compute I{a^,x,y, h) in (49), we apply the method of analytic stationary phase 
in the aa;-variable. 

Consider the auxiliary function 


T(aa;,ag,x,y) = ^{ax.,a^,x,y) - ^{a%,a^,x,y). 

Then, 9 q,^T(q;^, ag, x, y) = 0, 'k(a^, ag, x, y) = 0 with x, y) ~ (ag) and 

Im4'(aa,,a^,x,y) > 0 . 

Let C/(x,y) C C" be an open neighborhood of a^(x,y) on which the Morse Lemma 
holds (cf. |Sj2t Lemma 2.7] ), and set 

LR(x,y) = 17(x,y) nM”. (50) 

Further, let 5 G C'“(M"',M) be defined as 

(51) 

(««) 

for some di > 0 small. For each a^ G M”" let F = F(ag,x,y) : VM.{x,y) —)■ C"" be the 
complex contour given by F(a^,x,y) = Da^(zv^(^x,y)^ioix', x, y), where 

r{ax]a^,x,y) = Ox + i 6{a^) da^'if{ax,a^,x,y). 

We choose <5i small enough so that F(aa;; a^, x, y) C U{x,y). With this choice of 
contour, 

Im4'(F(aa;;a^,x,y),a^,x,y) = 6{a^)\da^'I>{ax,a^,x,y)\'^+0{6'^{a^)\da^'^{ax,a^,x,y)\‘^) 
for all ax G V. Since a%{x,y) is a non-degenerate critical point, 

ImT(aa;,a^,x,y) > G\ax - a%\^ 

for some C > 0 and all z G F(a^,x,y). Consider the boundary surface 5r(a^,x,y) : 

[0, (5(a^)] X 9Ci;(x,y) C" joining MR(x,y) and F(Q;^,x,y) given by S'r(a^,x,y) = 
Ute[o,5(a5)]5'r(L x, y), where 


Sr{t,ax;a^,x,y) := ax + it da^'i>{ax,a^,x,y). 
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Let Qr{ct^, x,y) C C” be the domain with boundary 

d^lr{a^,x,y) = T{a^,x,y) U VR(x,y) U Sr{a^,x,y). 



First, from 


one can write 


I{a^,x,y,h) = x, y;/i) + ^R"\yR(a5, y,/i), 

■JVr 

(52) 

where R^n\Y^{a^, x,y, h) denotes the integral over W^\V^{x,y). Then, by Stoke’s 
Theorem, 

I{a^, X, y, h) = J a^, x, y, h)dz 

+ Rr{a^,x,y,h) + RT^r.\Vj,{a^,x,y,h). 

(53) 

Here, the phase is T®^(z, a^, x, y) = ^'^{z,a^,x,y) — ^{a'^,a^,x,y), where denotes 
the holomorphic extension of the function d> defined in (47) in the a^-variable. In the 
first term on the RHS of (53), 

c^{z, a^, X, y, h) := a^{z, a^, y; h)b'^{z, a^, x, h)p^{z, y)p^{z, x), 

where p^ denotes the almost-analytic extension of p (see [Zwl Theorem 3.6]). The re¬ 
mainder i?r involves integration over Hr and 5r. An explicit description of Ry-, R^n\VK 
and all other subsequent remainder terms are given in §3.4.2 By choosing the con¬ 
tour deformation space Hr sufficiently small (after possibly rescaling the parameter 
t G [0,1])), it follows from the holomorphic Morse Lemma that there exist holomor¬ 
phic local coordinates w = {wi, ..., Wn) in a neighbourhood of Hp containing a^{x, y) 
such that 

T‘^(u;,ag,x,y) = i - - -(a^). 

Letting r(a^, x, y) be the image of r(a|, x, y) under the change of variables z w, 
one can write 

I{a^, X, y, h) = j g-^ <“s>c^(u;, a^, x, y, h)dw+{Rr+RRn\Vti)ioii^ y, h), 

(54) 

where c^{w, a^,x, y, h) := c^(w, a^,x, y, h) det (^(w, x, y)) . This choice of coordi¬ 
nates and the definition of r(a^, x, y) imply that 

Re[('u; - alf] > C\w - 
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for some C > 0 and all w G r(Q;^, x, y). 

By the Implicit Function Theorem, there exists an open snbset TQ{a^,x,y) C 
T{a^,x,y) and a neighborhhod Vf^{a^,x^y) C M” of a%{x,y) so that 

fo(a^,x,y) = {ax + iH{ax-,a^,x,y) : G V^^{a^,x,y)} 

where H{ - ■,a^,x,y) is an analytic function with \H{ax;a^,x,y)\ < A|a;a;| for some 
A < 1 independent of ax- Then, 

/(a^, X, y,/i) = [ e ^°'^'^c^{w,a^,x,y,h)dw 

Jto 

+ Rr{a^,x,y,h) + R^n\v^{a^,x,y,h) + R^-^f-^{a^,x,y,h), 

(55) 


where -Rr\fo comes from replacing the domain of the integral in (54) with r(x, y)\ro(x, y). 

As before, consider the boundary surface : [0,1] x —)• C” 

*^fo ^x) — Oix it Hi^axi a^, x, y) 

joining Tq with Fp^. Also, set np^(a^,x,y) C C"’ to be the domain whose boundary is 
To U Fp^ U S’p^. Then, another application of Stoke’s Theorem in ( |55[ ) gives 

I{a^, X, y, h) = f e“^ W)c(aa,, ag, x, y, h)dax 

K 

+ Rr{a^,x,y,h) + R^r.\v^{a^,x,y,h) + R^-^f-^{a^,x,y,h) + Rp^{a^,x,y,h), 

(56) 

where Rp^ involves integration over and S'p^. 

To finish the argument, we use that 

c^{z, ag, X, y; h) := a^{z, ag, y; h)W-{z, a^, x, h)f)^{z, y)p^{ax, x), 

and that (a^ • y^){z, a^,x, h) is holomorphic in z. Then, by standard asymptotics for 
Laplace integrals |Sj2 Theorem 2.1], 

r (aa;— \ 1 

/ e 2S ^°‘0[a-b){ax,a^,x,y,h)dax = Tz—r^aG{x,y,a^,h) + Qaa{x,y,a^,h), 

Jv- [Zirri) 

To 

(57) 

for ac G S^l ^ defined by 




aGix,y,a^,h) = ^ h^Wk{x,y,a^), 


k=0 


(a • 6)(a^,a|,x,y), 


where Cq is a positive constant. Here, 
Wk{x,y,a^) = ^ 


A 


k 


(58) 
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and 

\Qaaix,y,a^]h)\ < C (1 + i)2 e"^, 

for some C > 0. Note that, in particular, wo{a^,x,y) = (1 + 

Combining (55) with (57) gives 

I{a^,x,y,h) = + RQ{x,y,a^,h) 


with remainder 


{2'Khy 


Rq — Ry + R^n\Vu + '^r\fo + Rac- 


RaGix,y,a^,h) = 


Finally, it follows from (48), ( |46| ) and (60) that G{x,y,h) decomposes as 
G{x, y; h) = Ag{x, y, h) + Rq{x, y, h), 

where we have set 

X{x,y) 


and 


AG{x,y,h) = f 

{2Trhp J^n 


RG{x,y,h) = x{x,y) / R(j{x,y,a^,h)da^ 


for R^{x,y,a^,h) defined in (61). 

We now complete the proof of Proposition 
3.4.1. Leading term AG{x,y,h). Since (re, y) = exp^ 


pa:) 


we have 


(59) 


(60) 


(61) 


In (61), the additional remainder term i?i_p(a^, x, y, h) is given by 

f e"^^'W^<"e)(o • 6)(a,j,,a5,x,y)[l - y(aa:,x)y(aa,,y)]daa;, (62) 

K 

and 


(63) 

(64) 

(65) 

( 66 ) 


4>(a^(x,y),ag,x,y) = -2(exp„c (y), ag)«c + -dg(x, y)(a 5 )«c. 

Given p,q € M, consider the parallel transport operator (along the unique shortest 
geodesic from q to p) Tq-^p : T*M ^ T*M. This map is an isometry that satisfies 

Tq^p exp"^ (p) = - expp^ (y) and Tq^p = 7),%^. 

Changing variables ^ y := fac{x,y)^y{a^), where fo,c(^x,y)^y ■ 

^y Wx 


map induced by the choice of coordinates, and using that exp ^{a%) = i exp ^(x), we 


get from (65) that in local coordinates 

x{x-,y) 


denotes the 

2 


Ag(x, y, h) = ^ ^ x, y, h)dy, 

[2Trh)^ jRn 

V^(x,y,y) := (y“^ exp“^(x), y) + ^d^(x, y)(y“^/%). 


with 


(67) 
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and where after some abuse of notation we have set 


fda^ 


aG{x,y,ri,h) := a g{ x, y, Ty^(rj), h) det [-^id,x,y) 
for ac G defined in ( [5^ . In particular, since \Ty^ag{'n)\a^ = Ivly, we have 

ing term AG{x,y, h) in Proposition]^ 


_ In particular, since \Ty- 

^ det (^(r/,x,y)). This proves the identity (29) for the lead- 


3.4.2. Remainder term Rg{x, y, h). We proceed to prove statement (30) in Proposition 
In the notation of Theorem 


with 


G(x, y, h) = Ag{x, y, h) + Rg{x, y, h), 
Rg{x, y, h) = Rq{x, y, h) + R{x, y, h). 


Here, we recall that R{h) = —G{h)Rab{h){I + Rab{h)) ^ as defined in (45) and 

RG^x,y,h) = x{x,y) / RQ{a^,x,y,h)da^, 


where according to (61) 


Rq — Rr + R^-n\VR + -^r\fo -^fo + Ri-p- 


( 68 ) 


We now prove the exponential decay in h for each of the remainder terms comprising 
Rg{x, y, h). The exponential decay of d^dy Rg{x, y, h) is proved in the same way. 


Remainders R^n\y^ and R-p\f'^- The term (resp. R-p\f'^) is a result of shrinking 

the domain of integration M"' to V^{x,y) C M"' (resp. T to Tq C T). Namely, 

RR"\vA^,y,h) = x{x,y) [ [ ( 69 ) 

JR" JR^VMr 

To study the decay of R^n\y^, assume without loss of generality that I^(x, y) is a cube 
centered at a^(x, y) with side length 25o with do > 0 independent of x and y: 

PR(x,y) = {oj, G M" : - {al{x,y)Y'^^\ < 6o, k = l,...,n}. (70) 

Given ax G W^\Vu{x,y), we have that either dg{y,ax) > do/2 or dg{x,ax) > do/2. 
Consequently, since lm.^{ax,a^,x,y) = \{d?g{ax,x) + d?g{ax,y)){a^) 

d 2 

Im$(aa;,a^,x,y) > -^{a^)a^. 

Thus, there exists C'(do) > 0 with 

= 0 (e-^(^°)/"). 

The decay for is proved in the same way. 
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Remainders Rr and Rp^- The term R^ (resp. Rf^) is the result of an application of 
Stoke’s Theorem and consists of an integral over Sr and an integral over fir (resp. 
and More precisely, 

where 

RsYix,y]h) = x{x,y) [ [ e^'^^^^’°'^’^'y^c^{w,a^,x,y)dwda^, (71) 

JR" JSy 


and 

Rn^{x,y;h) = x{x,y) [ [ e^'^^^^'‘^^’'^’^^du,c^{w,a^,x,y)dwda^. (72) 


We first prove decay for Rsy- As before, for w G dVM,{x,y), either dg{w,x) > 6o/2 
or dg{w,y) > 5o/2. Also, by choosing Ji in (51) sufficiently small in terms of Jq; one 
can arrange that 

dg{Sr{a^,x,y), 514(x,y)) < 5q/A. 


By Taylor expanding ^'^{w,a^,x,y) at |Imt(;| = 0 we get that there exists (7 > 0 so 
that Im $‘i'(u;, x, y) > C{a^)new for all w G 5'r(a^, x, y). This gives the exponential 
decay of Rs^ ■ 

As for the remainder Rq^ {x, y, h), one uses the fact that the amplitude in the integral 
for Rfi^{x,y, h) contains the term dwC^{w, a^, x,y) and that 


dYoC^{w,a^,x,y) = {a ■ b)‘^{w,a^,x,y) ■ dw[p^{w,x)p^{w,y)]. 


We know that p^{Rew, x)p^{Rew,y) = 1 if both dg(Rew,x) < mi{M,g)/8 and 
dg(Rew,y) < mi{M, g)/8 hold. It follows that the integrand for Rn^{x,y, h) has 
its support contained in 


{tc G fir : dg(Reu;,a:) > inj(M, 5 r )/8 or dg{Rew,y) > g)/8}. 

The rest of the argument is the same as that for Rsj.{x, y, h). The analysis of the decay 
of Rf-^ is analogue to that of Rr so we omit it. 


Remainder i2i_p. The term Ri-p arises after removing the cut-off functions from the 
symbol so that the result is an analytic symbol and then one can apply analytic 


stationary phase for quadratic phase functions. It follows from (62) that 


Ri-p{x,y,h) = 

= f [ "2h"'‘ (“5)(a • b){aa;,a^,x,y)[l - p{ax,x)p{ax,y)]daxda^. 

JR" JVf, 

To 

(73) 


The integrand of Ri-p{x, y, h) is supported in the set of (a^,, x, y) G x M x M for 
which dg{ax,x) > m.]{M,g)/8 or dg{ax,y) > g)/8. Since the variable ax ranges 

over 1/p^, we deduce that Ri-p{x,y, h) = 0 unless dg{x,y) > Co for some Co > 0. The 
rest of the argument is the same as for i? 5 p(x, y, h). 


Remainder Ra^ ■ We recall from (63) that 
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Raaix,y,h) = 


i [ (ag) 


'Vf, 


/ / 

■/R" JR’ 




\^f 


[(a • b){ax, a^, x, y) - Ch{ax, a^,x, y)] duxda^ 

(74) 

Ch{ax,a^,x,y)daxda^, (75) 


ro 


where 


/ N 92,[o-6](a^,a^,x,y) 

Ch[ax,a^,x,y) := 2^ -- {ox - a^y. (76) 




7 ! 


3n n 
4 ’ 4 


The exponential decay of (74) follows from the fact that o G and b G 


3n n _ 9 

-* 4 ’ 4 


Indeed, the error term a ■ b — Ch, in the Laplace integral asymptotic (see (58) - (59)) 
satisfies the estimate 

|(a • 6 - Ch){ax,a^,x,y)\ < 

with Cl > 0. The exponential decay of (75) is obtained in the same way as for R^n^y^. 

Remainder R{x, y, h). Finally, we estimate i?(/i) = —G{h)Rab{h){I + Rab{h))~^. From 

(11, 

R{x,y,h)= / AG{x,u,h)Rabil +Rab)~^iu,y,h)du 
JM 


+ / R^{x,u,h)Rab{l + Rab) {u,y,h)du. 

JM 


(77) 


To deal with the second integral in (77) one simply uses the pointwise bound \Rq{x, u, h)\ 

(7(g-C//l) gg^ 


'M 


RMx,u,h)Rab{l +Rab) ^{u,y,h)du = 0{e ^1^). 


To estimate the first integral in ( [77| ), we note that |4 g'(x, tt,/i)| = 0(1), and use that 
the exponential decay of Rabix,y,h) in (42) to give R{x,y,h) = 0{e~^/^) uniformly 
for x,y € M. □ 
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